Irreversible decay of nonlocal entanglement via a "reservoir" of a single degree of 

freedom 
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Recently, it has been realized that nonlocal disentanglement may take a finite time as opposite to 
the asymptotic decay of local coherences. We find in this paper that a sudden irreversible death of 
entanglement takes place in a two atom optical Stern-Gerlach model. In particular, the one degree 
non dissipative environment here considered suddenly destroys the initial entanglement of any Bell's 
states \4> ± ) superposition. 



PACS numbers: 03.65.Yz, 03.65.Ud, 32.80.Lg, 42.50.Vk 



I. INTRODUCTION 

Superposition principle is a basic feature of Quantum 
Mechanics (QM). When applied to composite systems 
superposition leads to nonlocal coherences and quantum 
correlations (entanglement) which have been extensively 
studied in the last years, and are even expected of central 
interest for physics to come [lj|. 

As it is known, quantum coherences and entanglement 
are potential resources in various quantum information 
processes, even if they suffer any sort of environmental 
action which represents a serious obstacle towards these 
applications. On the other hand, it is just this fragility 
with respect to the environment, that is, the propensity 
to disperse information throughout inaccessible degrees 
of freedom that explicates the quantum-classical transi- 
tion @. 

An interesting point raised in [1, 0, H> S con- 
cerns the possible finite time disentanglement for bipar- 
tite systems, as opposite to the usual local decoherence 
asymptotic time. An experimental evidence of this pecu- 
liar trait of entanglement has been reported recently by 
Almeida et al. This issue has been analyzed Q in a 
simple and realistic model where two initially entangled 
two-level atoms separately interact with the multimodc 
vacuum noise of two distinct cavities. The authors find 
out that the nonlocal decoherence may take place sud- 
denly or at least as fast as the sum of the normal single 
atom decay rates. This sudden death of entanglement 
has been analyzed also for two Jaynes-Cummings (JC) 
atoms Also in this case the dynamics of the entan- 
glement between the atomic internal variables shows dif- 
ferent peculiarities for different initial states, with possi- 
ble sudden decays that are however followed by periodic 
revivals, due to recovery of information by the system, 
being the cavities lossless. 

In this paper we wish to investigate if the separate non 
dissipative interaction of each two-level atom with one 
only degree of freedom but of continuous variables, may 
play the role of the interaction with a reservoir leading to 



an irreversible disentanglement of the bipartite system. 

As well known, the optical Stern-Gerlach (OSG) model 
[10l | gives the opportunity of modifying the JC model by 
including the interaction between the internal and the 
external atomic dynamics via the e.m. mode of the cav- 
ity field, so actually providing the coupling of each qubit 
to one only degree of freedom. It has already been shown 
that this non dissipative coupling leads to a decoherence 
in the dynamics of a single atom [ll|, HH, EH and also af- 
fects the entanglement in the internal dynamics of two 
atoms that successively interact with the same cavity 
14j . As will shown in what follows, the same interaction 
gives rise to a disentanglement of the bipartite system 
which exhibits strong analogies with the disentanglement 
generated by the vacuum noise 0, Q ■ 



II. OPTICAL STERN-GERLACH MODEL FOR 
TWO ATOMS 



Let us consider two identical isolated two-level atoms, 
A and B, respectively crossing two distinct ideal cavities, 
a and b (in a direction orthogonal to the cavity axis), as 
in the OSG model. These two uncoupled subsystems are 
described by the total Hamiltonian 



H = H Aa + H 



Bb 



(1) 



where, in resonant conditions and in the rotating wave 
approximation for both atoms, 



H Aa = ^ + tk J j(a^a + Sf + \ 
2m \ 2 

+he sm(kx A ) ^S A + Sfa) (2a) 

H Bb = ?k + hu; (tfo + Sf + \ 
2m \ 2 

+hesm(kx B ) (b^S 3 + Sfb^j . (2b) 
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To simplify the notation we will sometimes use the in- 
dices A and B to indicate the two subsystems Aa and 
Bb, respectively. The cavity fields are, as usually, de- 
scribed by the bosonic operators a, a 1 * and b, b' , and the 
two cavity frequencies lu — ck, as well as the coupling 
constants e for the two subsystems, are assumed to be 
the same. The 1/2 spin operators S±, S{ account for 
the internal dynamics of j-th atom, while the conjugate 
variables Xj and pj describe the external motion of the 
same atom along the axis direction of the corresponding 
cavity (j =A,B). 

When the initial uncertainties Axj of the atomic wave 
packets are sufficiently small with respect to the wave- 
length A = 2ir/k of both modes a and b, the Hamiltonian 
([1]) assumes the form (see, for example, Ref. [HI for one 
atom case) 



H 



Pa 



2m 2m 



Pb 



+ ^L + hLU (N A + N B ) + hnA^ + m»fL» (3) 



tB -B 



where 



A* = 



N A 



1 



tfa + S* + -, (4) 



2VN A 

are constant of motion. Similarly for the subsystem B. 
Cl J x = 2ek\J~fFjXj is the (operatorial) Rabi frequency de- 
pending on the atomic position (with respect to a nodal 
point of the related cavity mode function) and on the 
excitation number of the j-th subsystem. Using previous 
results [HI for the one atom case, it is easy to show that 
the evolution operator may be written as 



U(t,0) 



-iHt/h 



exp 



2it 



ma^fi^XA 



x exp 



it 
2mh* 



Pa r ex P 



it 



a-N^PA 



x exp 



2 it 



ma^/tf x B > exp 



it 



■P% 



x exp 



it 



^nPxPb > < 



2mH 

M (t)(NA+N B ) 



(5) 



where 



ehk 



& N = ao^Nj, ao = , i9o(t) = u>t + ma t /(6h) (6) 

and it has been assumed that each atom enters the own 
cavity at the same time. 



own cavity and, consequently, of the coupling with its 
own translational motion. We will see that, differently 
from the usual JC model, in this case the elements of 
the reduced density matrices that describe the atomic 
internal dynamics show damped oscillations as a conse- 
quence of the OSG effect which causes a splitting of the 
atomic packets. As we will see later, this damping affects 
the entanglement decay whose rate is actually related to 
the phase space distance of the scattered packets. The 
section ends with a brief account of the results for the 
one atom OSG model, which are useful to explicate rela- 
tions and differences between the decay rates of nonlocal 
entanglement and local coherences. 



A. Two atoms 

Let us first suppose that the configuration of the entire 
system at time t = may be written as 



MQ)) = \<pa(P))\<pb(P))\x(o)) 



(7) 



where \<pa)i Wb) describe the translational dynamics 
along the cavity axes of atoms A and B, while 



| X (0)) = [cos 7 | + -) + sin 7 |-+)] |00) , (0 < 7 < - ) 

. . (8) 

shows that the two qubits are initially in a superposition 
of the Bell's states usually denoted |V' ± ), and the two 
cavities are in the vacuum state. The ket \-\ — ) indicates 
that the qubit A is in the upper state and the qubit B is 
in the lover state, and so on. Using the dressed states 



| X ±) = -L[| + 0).±|-1).] ( j = A,B (9) 



state (jHJ assumes the form 



| X (0)) = — {cos 7 [|xi> + |^>]|-0) B 

+ sin 7 [| X £> + |xi>] \-0) A }. (10) 

Applying the evolution operator ([5]) to the initial state ([7]) 
and using Eq. (JTUJ), one obtains the following expression 
for the state of the entire system at time t 



III. REDUCED DENSITY OPERATOR FOR THE 
QUBITS 

In this section we will analyze for the two different 
configurations considered in Ref.[§], the time evolution 
of the initial entanglement between the two qubits un- 
der the effects of the interaction of each qubit with its 



\m) = c ^e-^\Mt))i\<pm)\x + A ) 

+ \&))\Xa)]\-0)b 
+^e-«»W \ VA (jt)) \X B ) 

+ |0i(<)>|xi>]|-Q>A. 



(ii) 
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where we have set 



\<pf{t)) = ex P |T^"iao%| exp 



it 
Zmh 



\(fij(t)) = exp < -i 



2mH 



IVi(0)) 



(12) 
(13) 



and the following relations have been used 



(14) 



If the initial spatial distributions of the two atoms are 
gaussian functions of minimum uncertainty 



(Pj{xj,0) 



1 



1/2 



exp 



(xj - x jfi ) 2 
AAx 2 . 



(15) 



centered in Xjfi and with the same spread Axo, the x- 
representation of (|12p gives 



bf(x h t) 



exp ^ =p— maoXj 




x exp 



(16) 



where 



a;±(i) = Xjfi T a a t 2 /2, (3(t) = Ax 2 , + iHt/(2m). (17) 

As Eq. (fT3")) shows, \<fj(t)) describes the translational free 
motion of the j-th atom, while the ^-representation of 
\(j)f(t)) given by Eq. (fT6|) accounts for the well known 
splitting of the incoming wave packet in the OSG effect. 
Since we are interested in the entanglement of the two 
qubits, we will identify the internal atomic degrees as the 
system of interest, while the cavities and the external dy- 
namics will serve as environment. As a consequence, the 
density operator \ip(t)} {ip{t)\ must be traced on the cav- 
ity fields and on the atomic translation variables. Taking 
into account the normalization of kets (fT2|) and (fl3|) . and 
using relation © we obtain 



p A "{t) = Trtran*,field{\lp{t)) 

= oi | + -) (+-| + a 2 |-+) (-+| + a 3 

+^h+) (+-|+«4| — )(— I 



(18) 



where 



ai = 



«2 



«.3 



cos 2 7 [l + i?e(^(i)|^+(i)}] 
sin 2 7 [l + i?e(0 B (i)|0+(i))] 
cos 7 sin 7 «^(*)|^(*)) + 

x ((^(t)|¥> B (*)> + (^(*)IM*)>) 



a 4 = ---cos 2 7 .Re(4> A (t)|^+(t)) 
-isin 2 7 i? e (0 B (i)|^W) 



(19) 



By using the Eq.s (fTB"]) and (fH?)) it is pos sible to show 
that the following scalar products hold [l2l Il5| 

(0J(t)|#(t)) = eKp{-i%ot}e- r W 
<^(*)l0f (*)> = exp 
where 



ma?,t 3 fL nt 

i — =F « — — 

AH 2 



r(i) 



fe(i) 2 <5p(t) : 



8Ax 2 



8Ap 2 



(20a) 
-r(t)/4 (2Qb) 

(21) 



Apo is the initial uncertainty of the momentum distribu- 
tions for both atoms, and 

Sx(t) — — aot 2 , 5p(t) — —2maot, fij o = 2maoXj o/H. 

(22) 

In the following we will assume that xa,q = xb.q — x , 
which implies 

2m 



-a x . 



(23) 



By inserting this equation and Eq.s (|20p in Eq. p^)) we 
get 



ai = — cos 
2 



a-r = — sin 
2 



7 {l + cos(ft t)e~ r(t) j (24a) 

7 |l + cos(^ t)e" r(t) | (24b) 

a 3 = sin 7 cos 7 cos 2 (fV/2)e~ r( * )/2 (24c) 

a 4 = -{l-cos(Q *)e" r(t) }. (24d) 

It is useful to emphasize that T(t) is proportional to the 
square of the adimensional distance in phase space be- 
tween the average positions of scattered packets in the 
OSG effect, as defined in Ref. [l]| for the one atom case. 
The same T(t) can be also interpreted as the visibility of 
the Rabi oscillations [l2| in the ambit of the complemen- 
tarity relation as introduced by Englert [r| ■ 

Considering the ordered basis |++), |H — ), | — h), | ), 

the reduced density operator (|18p takes on the form 



/ \ 

ai B3 

S3 «2 

\ a 4 ) 



(25) 
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where (a\ + + CI4) = 1. 03 describes the coherence 
between the two atoms which assumes its maximum value 
at t = and, contrarily to what happens in the JC case, 
goes irreversibly to zero because of the coupling of each 
qubit with its own external dynamics. 

Now we consider the case in which the two qubits are 
initially in a superposition of the Bell's states denoted 
l^), and the two cavities are in the vacuum state, 

|x(0)) = [co S7 |++)+sin 7 |— }]|00). (26) 

With respect to the more convenient ordered basis |H — ), 

|++), I ), I — h), the reduced density matrix of the two 

qubits assumes the form 

( b 3 \ 

6 2 61 
b\ b 5 

\ 63 j 

where 

1 . 



AB 



(*) = 



(27) 



• sin 7 cos 7e 



-i0o(t) 



bi = 



6a = 



x((^(t)|^(t)> + (^ A (t)|0+(t)» 
x ((0i(t)b B (i)> + (4(t)|^(t)» 

1 + <^(t)|0j(t)) + i?e <^(t)|^(t)) 
+i? e (^(t)|0+(t))i? e (0 B (t)|0+(t))] 

1 + <£i(t)|^(*)) - i?e (<fe(*)l4(*)) 
-i2e<^(t)|^(t))flc<^(t)|^(t))] 

[l- Jte<^(t)|0i(*)) +Re(4>- B {t)\4> + B {t)) 

-Re(^(t)\^+(t))Re(^(t)\cl>+(t))] 

[I -Re {r A {t)\&)) ~ Re (4>g(t)\4>+(t)) 
+i?e(0 A (O|^W)i?e(0 B (O|^+W)]+sin 2 7 . 

(28) 

Using the results (|2T)|) we easily get 

61 = sin7cos7cos 2 (fi t/2)e _r(t)/2 

xexp{-i [tot- malt 3 /{3h)]} (29a) 

fe 2 = icos 2 7|l + cos(^o0e _r(t) } 2 (29b) 

63 = icos 2 7{l-cos 2 (f7ot)e- 2r(t) } (29c) 



cos 7 



cos 7 



cos 7 



cos 2 7 



65 = sin' 



7+ icos 2 7 {l -cos(£V)e~ r(t) j 2 (29d) 



where we have used 64 = 63 because of the condition 



B. One atom 

To analyze analogies and differences with the disentan- 
glement produced by different environments, it is useful 



to report in short on the dynamics of a single two-level 
atom in the OSG model. We assume that the initial state 
of whole system is 



IV>(o)> = Mo)) lx(o)> 



where 



|X(0)) = [cos 7 



sin 7 I 



10} 



(30) 



(31) 



and \tp) describes the atomic translational dynamics 
along the cavity axis. Using some results of Ref.pjl (see 
sec. IV) one obtains the time evolution of (f50|) . 



|V(t))=sin 7 Mt))|-0) 

+^e- WoW [k + (*)>|x + > + k-(*)>|x->] 



(32) 



Removing the environment by tracing on the field and 
on the atomic translation variables, one gets the reduced 
one qubit density matrix 



P(t) 



91 92 

92 93 



(33) 



where 



qi = _ cos 2 7 [1 + Re (4>-(t) I <£+(*)>] 



(34a) 



92 = -e- w «W cos7sin 7 [<^(t)|0+(t)> + (<p(t)\(f>- (t))] 



(34b) 



g 3 = -cos 2 7[l- J Re(0-(i)|0 + (i)>] +sin 2 7 . (34c) 



Both the irreversible evolution of the populations f (|34a|) 
and (|34c|) ) and the decoherence ()34b[) are to be ascribed 
to the increasing distance in the phase space between 
the scattered atomic packets. However, while the popu- 
lation decay depends on the square distance T(t) between 
the packets moving in opposite directions (through the 
scalar product {<jy^ l^)), the coherence experiences a four 
time slower decay since it depends on the square distance 
r(i)/4 between one of the moving packet and the stand- 
ing packet related to the ground state (scalar product 
(ifl^)). It is moreover to stress that, because of the 
non dissipative character of the decay, the qubit upper 
population does not go to zero, as in the spontaneous 
emission, but tends to half of the initial value. 



IV. ENTANGLEMENT EVOLUTION FOR THE 
TWO QUBITS 

A. Concurrence 

As it is generally convenient for these bipartite mixed 
states we will analyze the time evolution of the entan- 
glement between the internal variables of the two atoms, 
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Qot/zr Q t//r 



FIG. 1: Decay of concurrence, as given by Eq. (J35J) , when 
the two qubits are initially in the superposition (|HJ of Bell's 
states |V' ± )- The superposition parameter is 7 = 7r/4 (con- 
tinuous line), 7r/6 (dashed line), n/12 (dotted line), with 
xa,o = xb,o = xo — A/10, Axo — A/50. The other pa- 
rameters are A = 10~ 2 meter, e = 10 4 sec _1 , m = 10~ 2l 'kg. 



looking at the concurrence [17| . It is not difficult to show 
that the concurrence of the state (|25|) is 



C{p AB ) = max{0, (\\a 3 \ + Vai02l ~ Ikl - v^i^l)} 



= 2a 3 . 



(35) 



The entanglement with the non relevant systems causes 
a degradation of the nonlocal correlation between the 
qubits. According to previous results (see, for exam- 
ple 0, @), and as shown in Fig[TJ state (jHJ) gives rise 
to asymptotic disentanglement with the same rate of the 
decoherence. Because of the particular dynamics of the 
OSG model, in this case the disentanglement rate is one 
half time the population decay rate. Finally, we note 
that, for r(i) = 0, i.e. when translation dynamics of the 
center of mass of the two atoms is not taken into account, 
the concurrence (|35|) reduces to that derived in for the 
two JC atoms case. 

More interesting results are expected for the initial 
state (j2"rlj) . for which the concurrence is 



AB \ 



max {0, d(t)} 



where 



d(t) 



2(N -63). 
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(36) 



(37) 



The behavior of concurrence (|36|) is shown in Fig. 2. Also 
in this case, the coherence b\ between the two atoms 
goes to zero as a 3 in the previous case does. However, 
the concurrence experiences in this case a sudden death 
because of the probabilities 63 of finding the atoms in 
the two states |H — ) and | — h) which arc not correlated 
to the others and which tend to cos 2 (7)/4 for long time. 
Figure ||3J| makes evident the strong dependence of the 
entanglement life-time on the parameter of the initial two 



FIG. 2: Decay of concurrence, as given by Eq. (|36[) . when 
the two qubits are initially in the superposition (|26p of Bell's 
states \4> /■ The values of the parameters are as in figure 1. 



qubit state (126[) . similarly to the theoretical [J] and ex- 
perimental [8j results. However, in this two last papers 
it is outlined that only some initial entanglements may 
undergo to sudden death. On the contrary, due to the 
particular environmental action of the OSG model, there 
is no threshold in our case. In fact, the function d(t) of 
Eq. (|3T|) is negative if 63/I&1I > 1 and it turns out that 



& 3 /N>^sinh[r(t)/2]. 



(38) 



From this equation one can realize that for any value of 7 
there exists a finite time for which the concurrence goes 
to negative values and the two qubits become disentan- 
gled. In this respect, our model behaves similarly to the 
continuous variable two-atom model of Ref. [B| . 



B. Separability 

We will now inquire on the separability of the re- 
duced density matrices ([23]) and ([77)) adopting the Peres- 
Horodecky test 18, As it is known, for two qubits 
the non-negativity of the eigenvalues of the partial trans- 
posed matrix works as a necessary and sufficient condi- 
tion for the separability. Consider first the partial trans- 
posed of ([25]). 



AB 



(t) 



/ a* 3 \ 

ax 

a 2 

y 03 04 / 



(39) 



whose elements are given by Eq.s (|24[) . We easily find 
that three eigenvalues of (f3"9"| . 



Ai — ai, A2 — ci2, A3 — 



04 



4ai 



,1/2 



are non-negative, while the forth eigenvalue 



«4 

y 



1 



4a 



21 V2 



(40) 



(41) 



6 



. 4 
0.2 




-0 . 8 



1 2 3 4 5 6 
Q- t/n 

FIG. 3: Time behavior of the eigenvalue (|4ip of the matrix 
(|39l) for the one excitation initial state (jSJ . An enlarged form 
of this figure should show that A4 takes on the zero value 
only for isolated points. The values of the parameters are as 
in figure 1. 



0.4 



0.2 




-0 . 6 



-0 . 8 



1 2 3 4 5 6 

Q- t/n 

FIG. 4: Time behavior of the forth eigenvalue (|44[) of the 
matrix (|42[1 for the two excitations initial state (|26[1 , The 
values of the parameters are as in figure 1. 



takes on negative values for some time intervals (see fig- 
ure ((4])), in exact correspondence with the existence of a 
concurrence (see figure ©)■ 

V. CONCLUSIONS 

In this paper it has been analytically derived the time 
evolution of the initial entanglement of two atoms which 
separately interact with two distinct cavities. We have 
used the optical Stern-Gerlach model to include the ex- 
ternal translational dynamics of each atom. The pres- 
ence of entanglement has been ascertained trough the 
concurrence and the separability and we have found, as 
in Ref. [4|, that the two atom internal dynamics follows 
different patterns and becomes irreversibly separable in 
different times, depending on the different initial config- 
urations. As found in other contests [J, [fj, [8| which use 
different environments and where the disentanglement is 
due essentially to spontaneous emission, also in the OSG 
model the entanglement may undergo to an irreversible 
sudden death. In addition, the particular non dissipa- 
tive environment here considered suddenly destroys the 
initial entanglement of Bell's states \<p ) superposition 
(Eq. I26[) for any initial value of 7. This can in part be 
ascribed to the peculiar dynamics of the qubits in the 
OSG model in which the population 63 (matrix (|27[) ) of 
the non correlates states \-\ — ) and | — h) does not go to 
zero, but tend to cos 2 (7)/4 for long time. We can con- 
sequently conclude that nonlocal coherences manifest an 
amazing fragility even towards a non dissipative environ- 
ment made of a single degree of freedom. 



assumes only non-positive values, as figure §3§ shows. 

We wish to point out that for the initial state ([5]) the 
two qubits display EPR correlations [2(| for all times, 
but for isolated points (see Fig.([T])). Accordingly, one 
can see from the analytical expression of A4 (more accu- 
rately than from Fig. ([3])) that the bipartite system is 
not separable but for isolated points, and tends to the 
separability for long times. 

Finally, for the density matrix ([27)) relative to the ini- 
tial state (|2"rj)) the partial transposition gives 



AB 



(t) 



( h b\ \ 

6 2 

b 5 

V h b 3 J 



(42) 



whose elements are given in this case by Eq.s (|29]) . As in 
the previous case, three eigenvalues 

u 1 = b 2 , ^2 = 65, ^3 = ^3 + |6i| (43) 

are non-negative, while a forth eigenvalue 

vt = h-\bi\ (44) 
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